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Dark matter annihilation has been proposed to explain the TeV gamma rays observed from the 
Galactic Center. We study constraints on this hypothesis coming from the mass profile around 
the Galactic Center measured by observing stellar dynamics. We show that for several proposed 
WIMP candidates, the constraints on the dark matter density profile from measurements of mass by 
infrared observations are comparable to the constraints from the measurements of the TeV source 
extension. 
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A significant density of Dark Matter (DM) in the Universe has been observed on many length scales. The first 
evidence of the current dark matter problem came from the dynamics of the Coma cluster . Evidence of DM on a 
galactic scale came from rotation curves of galaxies which show that the orbital velocities of stars in galaxies do not 
follow the mass density derived from cataloging the luminous matter 0- This discrepancy can be resolved by adding 
a large amount of dark matter, a DM halo, that would not be included in a count of stars, gas and dust. More recent 
evidence of DM includes observations of the anisotropy in the cosmic microwave background, the luminosity-redshift 
relation for supernovae, and the theory of BigBangNucleosynthesis which measure a baryonic density of ili,h 2 = 0.022 
fNJ ' and a total matter density of VL m h 2 = 0.13 0, 0,1a, Ej. This implies that over 80% of the matter in the Universe is 
\ dark and non-baryonic. 

One class of candidates for non-baryonic DM is the weakly interacting massive particle, the WIMP. Theories such 
as supersymmetry, an extension of the usual space-time coordinates to include non-commuting coordinates, naturally 
include WIMP candidates. WIMPs are predicted to annihilate into other particles with energies similar to the original 
Oh' mass. The annihilation rate is a critical parameter in determining the relic density of these WIMPs and consequently 
one measure of whether they are a good candidate for a bulk of the dark matter. Photons will result from the 
annihilation, either directly, or through pion decay or acceleration of charged annihilation products. The annihilation 
\ could thus result in a "WIMP star" shining in gamma rays with energies near the particle mass. 

Recent advances in gamma-ray astronomy may allow the detection of DM annihilation. Ground based gamma-ray 
telescopes are currently sensitive to photons with energies above 100 GeV and have reached the sensitivity of a few 
percent of the Crab nebula flux. The most sensitive ground-based gamma-ray observatory currently in operation is 
/\ ' the HESS array in Namibia^- 

The Galactic Center (GC) has been proposed for observations of DM annihiliation HH[I3,[ll[ll[ll[ll[ll 
because it is close and might have a dense concentration of dark matter resulting in a strong signal of gamma 
rays. After tentative detections of a TeV gamma-ray flux from the GC by the VERITAS collaboration 

and the 

CANGAROO collaboration t ne HESS collaboration [l^ has initiated observations of the GC. HESS has reported 
a steady excess of TeV gamma-rays from the GC during two observational periods of 4.7 hours and 11.8 hours (at the 
6 sigma and 9 sigma levels respectively) . This excess of gamma rays is confined to a region of 3 arcminutes centered 
around Sagittarius A* , the dynamical center of the galaxy which is believed to host a supermassive black hole Q9| . 
The spectrum of this excess is a power law (4g ~ E~ a ) with a — 2.2 ± 0.2 in the energy range [0.2,8.8] TeV |20|. 
This spectrum is harder than typical gamma-ray sources such as plerions and active galactic nuclei. However, many 
of the newly discovered supernova remnants have similar spectra |2ll |22( . 

The GC gamma-ray flux may be produced by a variety of mechanisms p3l I24I I25II2I I . For example, the central 4 
million solar mass black hole could produce the gamma-ray flux by accelerating electrons in an extreme advection- 
dominated accretion flow, there is a high rate of supernovae near the GC and the shock fronts could accelerate particles 
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to TeV energies. Alternatively it could be due to DM annihilation. The GC has been suggested as a possible site of 
enhanced DM annihilation fl3L I2H l28l 12^. l30l | because it has a large stellar cusp and a million solar mass black hole 
in the center. The minimum radius to which any central dark matter density features extend is a key unknown in 
predicting the gamma-ray flux from WIMP annihilation. The interpretation of the photon flux from the GC is not 
settled and could lead at the very least to another example of an extreme particle accelerator, and possibly could shed 
light on the dark matter problem. 

In this paper we consider the DM interpretation of the GC gamma-ray emission and study if any more information 
on this hypothesis is contained in the dynamical measurements of the mass profile at the GC. Advancements in 
infrared astronomy are testing the small scale mass profile of the center of the Milky Way, down to tens of AU. With 
the WM. Keck 10m telescope, proper motions of stars have been monitored near the GC since 1996 pH l3ll l32|. 
Entire orbits have been or will soon be measured around the dynamical center of the galaxy. A strong gamma-ray 
signal from the GC implies a large amount of DM under the WIMP annihilation hypothesis. The change in mass 
enclosed in a sphere with radius d, the distance from the GC to the star, changes as the stars, on highly elliptical 
orbits, traverse any central spherically symmetric density enhancement of the dark matter. This could lead to an 
observable deviation of the orbit from a purely Keplerian orbit. Upcoming observations will provide direct constraints 
on the DM density profile in the center of the Milky Way (3j| and help us interpret the gamma-ray flux from the GC. 

We use the data on the stellar orbits around the GC published in [34|. More recent data, including the complete 
orbital shapes, may provide further constraints [3i| • We find that the gamma-ray flux from the GC is compatible 
with annihilation of a heavy, ~ 10 TeV, DM particle with a density profile consistent with the stellar orbits near the 
GC. Depending on the particle physics assumptions, the stellar orbits constraint is comparable but slightly stronger 
than the constraint on the source extension due to the angular resolution of HESS. Gamma-ray observations could 
have a very strong signature of WIMP annihilation due to the process \X ~ * 77: which would create a monochromatic 
line in the energy spectrum at the mass of the annihilating particle. Unfortunately, we find that for a TeV neutralino 
the flux of the monochromatic line is too weak to be seen with an energy resolution of 10%, the resolution of the 
atmospheric Cherenkov method. 

In the next Section we review the analysis used to connect the gamma-ray emission to the stellar dynamics consid- 
ered. We define expressions for the expected gamma-ray flux from WIMP annihilation with the purpose of clarifying 
its angular dependence and the units. Next we discuss the dark matter profiles we will use. We then study the limits 
imposed on dark matter at the GC by the astronomical mass measurements and the HESS angular profile. Finally, 
in Section III we present the conclusions of the analysis. 



II. ANALYSIS 



The flux of photons produced by DM annihilation depends on four factors: the annihilation products energy 
spectrum, the DM particle mass, its annihilation cross-section, and the density of the DM particles. The energy 
spectrum of the annihilation products, the annihilation cross-section, and the particle mass can be calculated once a 
particle model is specified. The density profile of a dark matter halo has long been a subject of much debate in the 
literature. Theoretical astrophysical considerations and numerical simulations have been used to suggest a family of 
DM halo shapes that could exist. The WIMP annihilation rate is proportional to the square of the particle density 
and many of the suggested DM halo shapes formally diverge when the emission rate is integrated along the line of 
sight through the center of the halo. The annihilation flux in the center of the DM halo will dominate the flux from 
these divergent halos. Astrophysically the density profile is expected to flatten at small radii where infalling objects 
can sweep out the centers of the halos through dynamical heating, although adiabatic accretion onto central baryonic 
density enhancements in the centers may create dark matter density enhancements |27], H^, |2^| . The annihilation 
rate is ultimately expected to limit the DM density [ilj. l30| . 

As a reference and to clarify the units of the quantities involved, we derive the expression of the photon flux from 
WIMP annihilation. Consider a small emitting volume dV at a distance I from a detector of collecting area dA 
(orthogonal to the line of sight.) This volume subtends a solid angle d£l as seen from the detector. Let dN e be the 
number of photons emitted during a time interval dt from the volume dV. Assuming the emission is isotropic, a 
fraction dA/(4ir£ 2 ) of the emitted photons is detected. Thus the number of detected photons in the same amount of 
time dt is 

dA 

dNn = dN e ^. (1) 
Specifically, for WIMP annihilation (xx ~^ anything — > 7) , the number of photons emitted is 

dN * = \ 3& "4- M dE dt dV - ( 2 ) 
2 da mi 



3 



Here (av) is the x~X annihilation cross-section times relative velocity, p is the WIMP mass density, m x is the WIMP 
mass, and dN 1 /dE is the number of photons in the energy interval [E, E + dE] produced per annihilation. The factor 
of 1/2 is there because 2 WIMPs are required per annihilation, p 2 {crv) dtdV jm 2 x is the number of WIMPs annihilating 
and dNj/dE is defined per annihilation. The photon flux from dV per unit energy at the detector then follows as 

rf£ dN D _J_^iP^ ( )r]v /q\ 
dE dAdtdE 8n£ 2 dE m 2 x ' ' ' 

When Eq. © is integrated along the line of sight, dV can conveniently be written in terms of the solid angle dil as 
dV = dlt 2 dVl. This leads to the usual formula for the flux per unit energy per unit solid angle, 

d$ 1 dAL (av) dJ 

(4) 



dEdVL 8tt dE m\ dQ 
where 

with the integral taken along the line of sight. We have written the solid angle d£l explicitly in dJ/dft to stress that 
its units are (mass density) 2 x (length) /(solid angle), as follows from Eq. and our derivation. This same quantity 
is denoted by J(tp) in the literature, e.g. yj . 

Eq. (0J can be integrated over a region 1Z of the sky to give 

d$ _ f d$ JO _ 1 dN 1 (av) f dJ 



dE J dEdfl 8tt dE ml J dfl 

n iz 



dft = - — — ? ^ / — dfl (6) 



When integrating over the whole source, Eq. JHJ gives a total flux of 



d$ _ 1 dN xx _ 



dE 8tt dE m\ 



(7) 



where 



J has units of (mass density) 2 x (length). Several units have been used in the literature. In particular Bergstrom, 
Ullio, and Buckley [12( used 8.5 kpc (0.3 GeV c~ 2 cm -3 ) 2 . For brevity, we introduce a Bergstrom-Ullio-Buckley Unit 
(BUBU) 

1 BUBU = 8.5 kpc (0.3 GeV c -2 cm" 3 ) 2 = 2.3605 x 10 21 GeV 2 c -4 cnT 5 = 0.530734 M% pc" 5 (9) 

Thence we will quote dJ/dfl in BUBU sr _1 and J in BUBU. These units were chosen so that a cored isothermal 
profile for the Milky Way halo would have dJ/dtt ~ 1 in the direction of the Galactic Center. 

For a source whose size R is small compared to its distance D, we can replace i in Eq. @ by the source distance 
and use cartesian coordinates centered at the source. We write the volume element dV — dx dy dz where z is along the 
line of sight and cc,y are transverse to the line of sight. To study the angular dependence of the signal, we integrate 
in z only and introduce the angles 6 X = x/D and 9 y = y/D. In terms of these, Eq. J3J gives 



d<£ 1 dN xx ^ (av) dJ 



dEd 2 9 8tt dE m 2 d 2 9 



(10) 



where 

d 2 9 

Integrating Eq. JHJl over the small source (R <C D) gives Eq. @ with 

1 r , 



' IJ = I p 2 dz. (I'll 



J =^2 I P dV - (12) 
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FIG. 1: Minimal supergravity models of the annihilation flux fit to the 2004 HESS data. The grey band consists of several 
spectra generated with mSUGRA neutralinos of mass ~ 11 TeV. 



A. Particle Model Examples 



Particle physics enters the gamma-ray flux through the combination 

dN 7 (av) 
dE ml 



(13) 



in Eq. J3J). We can estimate values for 4^, (crv), and the particle mass m x in examples of WIMPs. Once these values 
are given in a specific model, the resulting normalization required to fit the spectrum to the HESS flux gives a value 
for J. Varying the model parameters results in a band of J values. 

We give here three examples of WIMPs: the lightest neutralino in minimal supergravity (mSUGRA), the lightest 
neutralino in a generic minimal supersymmetric standard model (MSSM), and a Kaluza-Klein (KK) dark matter 
particle |26j| . 

To explore mSUGRA models we used the program DarkSUSY 35] to find model parameters consistent with particle 
accelerator and direct search bounds. The spectrum of gamma rays extends up to ~ 9 TeV and any WIMP annihilation 
that would explain the observation would require a particle with a mass above 10 TeV. In mSUGRA excessive thermal 
relic densities are predicted for most neutralinos with such a high mass. However, changing the cosmological model 
may alleviate this difficulty pjfl l37|. so we proceed without imposing the usual relic density constraint. We fit the 
normalization of the spectra to the HESS data. The results are shown in Fig. QJ. The two physical processes included 
in this spectrum are secondary pion decay and direct annihilation into photons. The spectral line due to direct photon 
production is not observable in the spectrum after it has been convolved with the HESS energy resolution of ~ 15%. 
Other processes, especially the acceleration of charged secondaries, could be reasonably expected to alter the spectrum 
|38| . This could provide other signatures of the annihilation which could be an important check on the DM annihilation 
interpretation of the HESS flux. We find that there is a family of mSUGRA models that produce a neutralino with 
a mass of 10 to 11 TeV consistent with current constraints and have a decent agreement with the HESS spectrum, 
with a x 2 °f ~ 1-2. These models require the J parameter to be in the range [300,3000] BUBU to explain the flux 
observed by HESS. 
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Lower values of J may be obtained once the parameter space is relaxed beyond mSUGRA. The difficulty in finding 
mSUGRA models that fit the HESS data lies in the excessive thermal relic densities predicted for neutralinos with 
the required mass, ~ 10 TeV or higher. Profumo 39} has suggested that resonant annihilation of neutralinos through 
the A boson in the early Universe, which can occur for to a — 2to x , can lower the relic density for particles around 10 
TeV. In this case the value of J can be as low as ~ 1600 BUBU (see his Fig. (7b) obtained in an anomaly mediated 
supersymmetry breaking model) or even ~ 3 BUBU (see his Fig. 8b, for a generic MSSM model). 

A third example of WIMPs that fit the HESS data are Kaluza-Klein particles. Ref. [2(| finds that the spectrum 
for a model of KK DM requires a J value of ~ 130 BUBU in order to be responsible for the flux recorded by HESS. 



B. DM density profile 

A dark matter density profile that would explain the TeV gamma-ray flux from the Galactic Center with the 
particles discussed in the last section will need to have a higher than expected density. The HESS GC source does 
not extend beyond ~ 0.1° or ~ 10 pc, covering a solid angle fi ~ 10~ 5 sr. The required J value from section III Al 
ranges from 3 - 3000 BUBU or a dJ/dVl ~ 10 5 - 10 8 BUBU sr" 1 for 0.1° or 9 - 10~ 3 rad. 

The contribution from the extended DM halo along the line of sight to the GC can be estimated from Eq. JSJ as 
dJ/dtt ~ p 2 D. For a canonical isothermal halo p ~ 2 pi oca i ~ 0.6 GeV cm -3 and dJ/dQ ~ 4 BUBU sr _1 from the 
DM column through the extended DM halo. This is 5 to 8 orders of magnitude smaller than required, so higher DM 
densities are needed to produce the gamma-ray flux by annihilation of our candidate DM particles. 

A Navarro- Frenk- White profile (NFW) |4jj is denser at the center. For an NFW profile, the average value of dJ /dVL 
within 0.1° is, from Eq. (|16J) below, 

dJ 1 f dJ 2itp 2 r 2 o , . , 

1F; dn = r^k ~ 3 x 10 BUBU sr ■ ( 14 ) 



dQ n J n dfl DO 

This is still 2 to 5 orders of magnitude too small to explain the observed gamma ray flux for most of the dark matter 
particles we consider. We conclude that, if the HESS signal is due to DM annihilation, the extended halo contributes 
no more than a few percent of the gamma-ray flux and a strong enhancement in the density must exist within 10 pc 
of the center of the galaxy. 

The dark matter density profile within 10 pc of the Galactic Center is not known in detail and mechanisms for such 
density enhancements have been proposed. For example, such enhancement could be explained by extreme clumping 
of the dark matter [ill l42l l43l |44| which would have implications on models of structure formation, by steeper density 
profiles, p oc r~ a with a > 3/2, have been su gges ted k5j but are disfavored, or by a strong dark matter concentration 
at the Galactic Center (a spike 0, [53, 128L l29t H(],E|). To include the latter two possibilities, we split the dark 
matter profile into an inner and an outer part at a transition radius Rj. 

As an example of the outer profile we use the NFW profile 



Pnfw 



1 + 



(15) 



r s is a scale radius and p s is twice the density at r s . We will take typical values [12] of r s = 25 kpc and p s = 
po(D/r s )[l +JD/r s ) 2 ] with a local density po — 0.3 GeV cm~ 3 . We take the distance to the Galactic Center to be 
D = 8 kpc |47|. For this profile we compute 



dJ NFW 2 J 7r — 9 3 + 2y 2 



P 2 s r 



dn rsh ) y 2(l+y2) 



arctan 



2(1 + 2T 2 )(1 + 7/ 2 ) 



(16) 



where 9 is the angle between the line of sight and the GC, x — D/r s , y = xs'm9, and z = xcos9. 

Notice that Eq. Ijl6|l diverges in the direction of the GC (9 = 0, y = 0) as irp 2 s r 2 /D9. To remove the inner part of 
the NFW profile, we add an inner cutoff at R\ by replacing the term (it — 9)/y in Eq. Ijl6|l with 

Hy, %) + HV, h) {l + y ^ (l + xl) arctan I i , 17) 



where x c = Ri/r s , z c = \Jx 2 - y 2 , b c = (zz c + y 2 )/{z - z c ), and 

1 p(l hl._ 

^TW 2 V 2 ' 2 ' 2 ! y 2 + aV v v » 



T{y,a) = -^= 2 F X (\,\^-^—^ =iarctan(^) (IM 
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The form with the hypergeometric function is used to avoid division by zero at y = (9 = 0). 

In the inner (r < 10 pc) part of the density profile we use a simple functional form to model a central density 
enhancement. We assume that a DM mass M\ is contained within a sphere of radius R\, and that its density profile 
is spherically symmetric and decreases with a power a of the radius. The inner profile we use is 



' 3 - a Mi 
Pi(r) = { if V^i 

,0, 



Rc < r < Ri, 
otherwise. 



(19) 



This inner profile could be a steep profile, or a spike around the central black hole. For this inner profile, we find 

(3 -a) 2 Mf 2 1 



■h 



3-2a 



4tt RfD 2 3 - 2a 

Here R c is an inner cutoff radius discussed in the next paragraph. For the angular profile we compute 



(20) 
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Y{a) 



(21) 



<9<9 h 



and zero for 6 > 9j. Here we defined 9\ = Ri/D, 9 C = R c /D, and 

F(a, x 



2 Fx(a,l;a + -;x ) 



(22) 



where 2 -Pi is the hypergeometric function. Notice that for a = 3/2 the factor in front of the square bracket is 
[27rln(6'i/0 c ] _1 and that for a = 1/2 the square bracket is \n(8i/9 c ). 

An inner cutoff at R c is introduced to avoid the divergence that occurs in dJi/dft when a > 3/2. This inner cutoff 
is left as a free parameter, because this part of the halo is even more unknown than the rest. Physically an inner 
cutoff would naturally be present. Either the capture radius of the black hole, or some effective radius at which, e.g., 
the DM density is depleted by annihilation during the history of the Milky Way. In the latter case the maximum 
sustainable density is usually taken as 



P < Pr, 



m 
(av)t 



(23) 



with the time t taken as the age of the Milky Way. In the case we are considering, we have a measurement of the 
flux and of the particle mass from the extent of the spectrum. For example, integrating Eq. Q in energy above the 
threshold and inserting Eq. (|19f) with a = 0, Eq. I|23|l implies 



Mi > M c = 



8irD 2 mt<I> 



(24) 



where $ is the total photon flux above threshold and iV 7 is the number of photons produced above threshold in each 
annihilation. Thus the maximum density p m ax corresponds to a lower limit on the mass of an inner feature of the halo 
that could explain the gamma-ray observation. If the mass is too small, then the cross-section and density required 
to mantain the same flux are so large that the feature would have annihilated by now. This can be generalized to all 
a values by finding the R c for which p(R c ) < pmax- They are any R c greater than the solution for R q in the equation 



r = r "[ 1 + Wj^rJ) ■ (25) 

Another scale in this problem is the capture radius of a 3 million solar mass black hole, expected to be in the center 
of all of these profiles. We find that the capture radius, ~ 10 -7 pc, is greater than all R q . 
Thus, as a physically motivated number, we take the range of cutoff radii to be 

10~ 7 pc < R c < R h (26) 
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C. Limits from the HESS angular profile 



The angular distribution of photons in the HESS detector carries information on the source profile. Here we 
investigate the constraint on the source profile due to these data. 

The HESS analysis [2(j assumes a gaussian source profile, and gives a limit on the source angular size equal to 
< 3'. To determine the limit on our power-law sphere in Eq. (|f 9(1 . we compare the emission profile, Eq. (|21|) .with the 
angular distribution of detected photons. Fig. (2) in [2(J gives the photon counts Ci and their errors Si in each 9 2 
bin. Here 9i is the angle between the photon direction and the center of the excess. The center of the excess agrees to 
the position of the GC to well within the systematic errors in the pointing of the HESS array. The intrinsic angular 
profile 



dJi 



dJ Q 

In 



is convolved with the point spread function (psf ) of HESS as given in |2f 

1 a 2 



fpsf(9) - fo 



T i H e 

8.7 



(27) 



(28) 



Here fo was chosen so the psf has unit area and the widths of the gaussians are a\ — 0.052° and 02 = 0.136°. We 
rewrite this as a linear combination of two normalized gaussians, 



re ' 



C2 T 

e " 



27rcr! 



(29) 



with ci = 8.7cr 2 /(8.7cr 2 + a 2 ) and C2 = <j 2 /(&.7<j 2 + cr 2 ). The source profile convolved with a normalized gaussian is 

i(f) 2 



and j convolved with the entire psf is 



a- 2 e 



3 P sf{6) = cii(6»,c7i) +c 2 j(0,a 2 )- 



(30) 



(31) 



The photon counts C (9) as a function of angle 9 from the GC are proportional to the convolution of dJ/dft with the 
psf 



C(9) = A 



dJ 



psf 



The proportionality constant is given by 



.4 



(32) 



(33) 



where £ is the exposure, N 7 is the total number of photons above the experimental threshold emitted per annihilation, 
and A9 2 is the aperture of the observation. We have estimated the HESS exposure as the ratio of the total counts 
assuming a point source and the integral flux of 1.82 x 10" 6 m~ 2 s^ 1 above threshold, both as reported by HESS in 
pfij . We estimate an exposure £~3x 10 13 cm 2 s. Furthermore, 



The best fit for the normalization factor A is 



.4 



dJ\ 

dn) 



E 



psf 



(34) 



(35) 



with the x 2 given by 



E 



Ci — AJjp S f{ 



(36) 
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FIG. 2: The constraint from the HESS angular size. The profiles are generally required to be steep, a > 1.5, or small, R\ < 1 
pc. 

Here to find our intervals we perform a bayesian analysis. We take the likelihood as proportional to e _x ' and 
define our confidence intervals in (a,i?i) as the corresponding quantiles of the posterior distribution. We take the 
prior distribution as flat in logi?i and a, and zero outside the range shown in Fig. |[2J. 

The intermediate results of this piece of the analysis are shown in Fig. J3J). The 68%, 90%, and 99% regions 
in the Rj,a are shown. At the 90% confidence level the HESS data confines the source diameter to < 1 pc for a 
uniform sphere. For power law density profiles with index a > 1.5 the constraint on the source size starts to weaken 
considerably; these profiles could be modeled as a smaller source with a harder power law index. Finally, we include 
a fit of two profiles in Fig. 10) • The first, more shallow profile is an NFW profile alone. Evidently the flux rules out 
an NFW only DM profile, and so does the angular profile of the gamma-ray source. The steeper profile is a fit of an 
NFW profile with a spike. A spike can simulate a point source effectively, and the extended NFW piece makes the 
profile fit slightly better than a point source. 

As an illustration of the kind of angular gamma-ray profiles compatible with the HESS data in models with a DM 
spike we show such a profile for an NFW model. We fix the halo NFW model with the parameters give in section 
111 131 A pure NFW model does not fit the data well at small angles from the GC. However, motivated by we a dd 
a spike with a radius of 1 pc and a slope of a = 1.9, which are typical values after inclusion of the effects of stars 
and annihilation [2£j. The \ 2 P er degree of freedom for the fit is 1.3 which is the same (A^ 2 ~ 10~ 4 ) as in the fit of 
a point-like source. Since the density of the dark matter is fixed for this profile is fixed to the local density and the 
particle mass is bounded by the spectrum, the normalization of the fit to the flux gives the value of the cross-section, 
N 1 {av) = 3.3 x 10~ 28 cm 3 s _1 . This is a reasonable value for the WIMP annihilation cross section in particle physics 
models. 



D. Limits from stellar dynamics mass measurements 

Measurements of the amount of mass M{< r) contained within a distance r from the GC are continuously improved 
as more precise data are collected. Here we use the compilation of enclosed mass measurements in |3 ll . From their 
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FIG. 3: An NFW+spike fit to the 2004 HESS gamma-ray data on the Galactic Center. For this fit we used a NFW profile 
outside 1 pc and a p ~ r -1 ' 9 profile for r < 1 pc. An inner cutoff of 10 -5 pc was used to approximate the inner core. This is a 
typical profile of a spike after dynamical heating and annihilation effects are included. 



paper, we extract a table of the mass Mj contained within radius fj, together with its quoted error G{. 

To these data we fit a mass profile with three components: the central black hole, the central stellar cluster, and 
the dark matter sphere described above. 

M(< r) = M BH + M*(< r) + M r (< r). (37) 

Mbh is the mass of the central black hole Sagittarius A*. For the central stellar cluster we use the empirical mass 
profile M*(< r) obtained from data in p8|. 




M*{<r) = l ) Jl *(i.o (38) 



with M* = 0.88 x 10 6 M Q and i?» = 0.3878 pc. We model the dark matter with the density profile described in 
Eq (|19|l . which corresponds to a DM mass profile 

r 



Mi«r) = <! Ml UJ ' (39) 
.Mi, r>R 1 . 

We use the likelihood function to find constraints on the DM density profile, using a bayesian analysis similar to 
section Til CI Assuming the errors quoted in |34| are gaussian, the likelihood £ is given by 



21n£ = V (Mi ~ M(< n))2 - V (Mi ~ MBH ~ M * (< n) - Mlfl)2 (40) 
11 ^ a? ^ a? 
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with 

/<=((S) ' n - Rl > (41) 

U, n>Ri. 

In order to obtain a constraint on the parameters M\ and i?i at a fixed value of a, we first marginalize over Mbh- 
Since In C is quadratic in Mbh , we need only replace Mbh in Eq. (|4U|I with the value M BH obtained by maximizing 
the likelihood. This is given by 



M BH = fgi^i, (42) 



with 



^i ^ Mj-MJ^n) ^ MM t ~M4<n)} ^ ft 

:r 1= ^-2, x 2 = }_^ -2 , Z3=2^ = £4 = 2^ ^ » x 5 = 1^^2- ( 43 ) 

As our prior, R\ is restricted to the range [0.0004,10] pc and distributed logarithmically, and a is kept at a few 
fixed values (0,1,2). By integrating our posterior probability distribution we derive a 1 sigma upper limit and a 90% 
bayesian interval in the M\, R\ parameter space. For R\ smaller than the innermost data point (0.0004 pc), there is 
a degeneracy between Mbh and M\. We break this degeneracy by imposing M\ less than the upper bound on the 
black hole mass reported in 34] (3.6±0.4 x 10 6 M©). This is equivalent to assuming all the mass within the innermost 
orbit could be DM. 



III. RESULTS 

From the particle examples in Section III Al we find that a range of J = [300, 3000] is needed to explain the flux 
of gamma-rays from the GC as DM annihilation products. With resonant annihilations, J can be as low as ~ 1. 
Furthermore we conclude that the DM annihilation line will be unobservable with an energy resolution of 10%. We 
find that the spectrum of gamma rays from the GC is compatible with the decay of pions produced in the annhilation. 
More complex models of the radiation, such as [3^] where bremsstrahlung of W products has an appreciable effect, 
the spectrum may be similar to a power law and other spectral features, such as a hardening of the spectrum near 
the WIMP mass, may be observable. 

The requirement that the central feature of dark matter not annihilate in the lifetime of the Universe gives a lower 
limit on the mass. For example, a central feature with a = and an upper limit on the density of p = 10 15 Mq pc -3 
and the requirement that J = 1000 BUBU gives a lower limit of ~ 3 x 10 -4 Mq. For a limit density of p = 10 12 Mq 
pc -3 the mass of annihilating dark matter must be greater than ~ 1 Mq to be stable for 10 10 years. These limits are 
below the lower edge of our results plots. 

The stellar dynamics limit extended mass distributions to ~ 10% of the black hole mass for R\ = (10~ 3 , 1). The 
angular size bounds are complementary excluding regions above a radius that depends on the assumed a for the 
distribution, as seen in Fig. (J2J. 

The results of the analysis are compiled in Figs. l0J[Slini and[7|). In Figs. JU[SJ[5Il we plot both the stellar dynamics 
bound and the angular size bound in the Mi — Ri plane for three values of a. The expected range from the particle 
physics are shown as shaded regions. These regions correspond to either a value of J that could produce the observed 
flux, or equivalently to a value for N(av). 

Comparing these models to our stellar dynamics bounds, we see that for a = (Fig. the source size is restricted 
to < 20 pc in mSUGRA and KK models. For larger cross-sections with resonant annihilation, the source size is 
unbounded by the stellar dynamics. The constraint from the HESS source profile limits the source size to < 1 
pc (vertical line), so it is similar to the stellar dynamics constraint in mSUGRA and KK models, but is stronger 
for resonant-annihilation models. However, for some of the mSUGRA models we considered the stellar dynamics 
constrains were stronger restricting the source size to < 0.3 pc. 

For profiles with shallow cusps (a = 1; Fig. O, the source size constraint on WIMP models from stellar dynamics 
is similar to the a = case. No bounds for resonant-annihilation models, but still < 20 pc for mSUGRA and KK 
models. The HESS constraint from the angular size of the gamma-ray excess is still ~ 1 pc and so conclusions similar 
to those with a — apply in this case. 

For profiles with steep cusps (a = 2; Fig.[f)|) stellar dynamics bounds out to 10 pc do not provide a constraint on 
the WIMP models we examined. The constraint from the HESS angular profile is also much weaker here. We show 
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FIG. 4: Stellar orbit constraints on the mass and size of a dark matter concentration at the Galactic center. The black lines 
are the 1 sigma upper limit and 90% confidence interval for the mass of dark matter spheres with power law density profile 
index a — 0. The dashed lines show the mass corresponding to the J values indicated (in BUBU). The 1 sigma and 90% bounds 
from the angular profile are shown as the vertical gray lines. The grey bands show the typical values of J required to produce 
the HESS flux in WIMP models: KK (dark grey), mSUGRA (medium grey), and MSSM (light grey) 



two plots here to illustrate the effect of the cutoff radius which only comes into play for these steep profiles. We show 
two cutoff radii of 1CT 4 and 10~ 6 pc. 

The constraints from stellar orbits and the HESS angular distribution are summarized for comparison in figure 
13 The solid line represents the 90% confidence region based on the HESS data alone. The dotted lines show the 
constraint coming from stellar dynamics. Various values of J are plotted so that these constraints can be compared to 
particle physics models. The values of J required by the mSUGRA neutralinos, MSSM neutralinos, and Kaluza-Klein 
particles we examined are plotted as the medium grey band. Both the stellar dynamics and the gamma-ray angular 
profile point to a DM source that is either small or steep. 

IV. CONCLUSIONS 

To summarize, a very high energy gamma-ray flux from the center of the Milky Way was significantly detected by 
the HESS collaboration during 2003-2004. A possible explanation of the very high energy radiation from the Galactic 
Center is WIMP annihilation. The intensity of the annihilation flux is a function of the density profile of dark matter 
in the galactic center. The angular distribution of detected gamma rays limits the size of the emission region. Data on 
on the proper motions of stars and star counts around the galactic center constrain the size and the mass of the dark 
matter at the GC. We have shown that the density needed to produce the observed flux from WIMP annihilation is 
consistent with observational constraints on the mass profile of the GC. For the stellar orbit data and the star counts, 
we used the infrared data in [34| and |4g. We found that these astronomical constraints on the source profile are 
comparable to and slightly stronger than the constraint from the angular distribution of photons measured by HESS. 

There are several ways in which WIMP annihilation as the origin of the HESS flux could be confirmed or made 
implausible. As is clear from fig.0 a slight improvement in either the gamma-ray angular resolution or the constraints 
from stellar orbits may reveal the presence of an extended dark matter annihilation region at the Galactic Center. An 
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FIG. 5: Same as Fig. |g) but for dark matter spheres with power law density profiles with a = 1. 



extended emission out to large angles would be a possible indication of WIMP annihilation. An extended gamma-ray 
excess with the same spectrum and position of the GC flux has recently been reported by HESS . A spectral cutoff 
at energies higher than the particle mass is another requirement of the DM hypothesis. The cutoff may be preceded 
by a gamma-ray line at the particle mass, but this spectral line does not appear to be observable with atmospheric 
Cherenkov telescopes in the particle models we examined due to the insufficient energy resolution. Absence of 
variability is another feature of WIMP annihilation, thus variability of the source would be difficult to reconcile with 
the DM interpretation of the GC TeV flux. Finally, since the dark matter permeates our Universe, if the same 
radiation was found in the centers of other mass concentrations, population studies may be possible that could help 
confirm or deny the annihilation nature of this radiation |Hcj . 

A small spike on an NFW profile could explain the large gamma-ray flux which is not expected from cored or cusped 
halos. Astrophysically small spikes in the DM halos are not favored, but not ruled out either. The infrared data of 
proper motions in the GC show about three million solar masses confined to a space of 90 AU. The compression of 
this baryonic matter may adiabatically compress the dark matter and lead to such a spike in the profile [13(. Any 
merger events with larger stellar sized objects should dynamically heat the DM spike reducing its density. 

Further observations of the Galactic Center in gamma rays are ongoing. There are hints that the TeV radiation 
from the Galactic Ridge is connected to the Galactic Center point source. The TeV flux from the GC seems to be 
constant in time and a cutoff in the spectrum (now reported to have a spectrum with a — 2.4) has not been found up 
to energies of ~ 6 TeV [5l|. so the models considered here are still viable. The nature of this non-thermal radiation 
source in the center of the Milky Way is still unknown and undergoing active study and observations. 
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FIG. 6: Same as Fig. J3) but for inner DM profiles with a — 2. For this steep density profile, which "looks" like a smaller 
sphere, the size constraint on the DM feature is not as strong, constraining values of J ~ 10 6 BUBU or higher. To avoid a 
divergence in J and because these profiles are expected to flatten at some inner radius, these profiles were computed with a 
minimum cutoff radius R c = 10 -4 pc (top) and R c = 10 -6 pc (bottom). 
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FIG. 7: Constraints on the J parameter due to stellar orbit data 34. 48] and the HESS source profile. The dashed lines show 
the crossing of specific J values with the 90% confidence level from stellar orbit measurements as shown in Figs. 14151 and [(J 
for specific values of a. The 90% constraint on the source profile from the HESS data alone is shown as the solid line. The 
grey band is the expected range of J for mSUGRA neutralinos that produce the GC TeV flux. 



models with heavy neutralinos. 



[1] F. Zwicky, Helv.Phys.Acta., 6, 110 (1993). 
[2] H.W. Babcock, 1939LicOB, 19, 41B (1939). 
[3] C.L. Bennett et al, ApJS, 148, 1 (2003). 

[4] R.H. Cyburt, B.D. Fields, K.A. Olive, Phys.Lett., B567, 227-234 (2003). 
[5] J.L. Tonry et al., Astrophys.J., 594, 1-24 (2003). 
[6] D.N. Spergel et al., astro-ph/0603449 ;2006) 

[7] W. Benbow for the H.E.S.S. Collaboration, 2nd Int. Symp. on High Energy Gamma Ray Astronomy, Heidelberg, APS 
Conf. Proc. 745, p. 468 (2004). 

[8] J. Silk and H. Bloemen, ApJ, 313L, 47S (1987). 

[9] F.W. Stecker and A.J. Tylka, ApJ, 343, 169S (1989). 
[10] R.A. Flores and J.R. Primack, ApJ, 427L, IF (1994). 
[11] H.A. Mayer-Hasselwander et al. A&A, 335, 161M (1998). 
[12] L. Bergstrom, P. Ullio, and J. Buckley, Astropart.Phys. 9 (1998) 137-162. 
[13] P. Gondolo and J. Silk, Phys.Rev.Lett. 83 (1999) 1719-1722. 
[14] A. Cesarini et al., APh, 21, 267C (2004). 

[15] N.W. Evans, F. Ferrer, and S. Sarkar, Phys.Rev. D69 (2004) 123501. 
[16] K. Kosack et al., ApJ, 608, (2004) L97. 
[17] K. Tsuchiya et al., ApJ, 606, (2004) L115. 

[18] W. Hofmann, Proc. 28th ICRC, Tsukuba, Univ. Academy Press, Tokyo, (2003) p.2811. 

[19] A.M. Ghez et al., ApJ, 620, (2005) 744. 

[20] F.A. Aharonian et al., A&A 425 (2004) L13-L17. 

[21] F.A. Aharonian et al., Nature 432 (2004) 75. 



15 



[22] F.A. Aharonian et al., A&A 432 (2005) L25-L29. 

[23] A. Atoyan, C. Dermer, ApJ, 617 (2004) L123-L126. 

[24] F.A. Aharonian et al, Astrophys.J. 619 (2005) 306-313. 

[25] D. Horns, Phys.Lett. B607 (2005) 225-232; Erratum-ibid. B611 (2005) 297. 

[26] L. Bergstrom, T. Bringmann, M. Eriksson, M. Gustafsson Phys.Rev.Lett. 94 (2005) 131301. 

[27] P. Ullio, H.S. Zhao, M. Kamionkowski, Phys. Rev. D, 64, (2001) 043504. 

[28] D. Merritt, M. Milosavljevic, L. Verde, and R. Jimenez. Phys. Rev. Lett., 88, (2002) 191301. 

[29] O.Y. Gnedin, J.R. Primack, Phys.Rev.Lett. 93 (2004) 061302. 

[30] G. Bertone, D. Merritt, Phys.Rev.D 72, (2005) 103502. 

[31] A.M. Ghez, B.L. Klein, M. Morris, and E.E. Becklin, ApJ, 509, (1998) 678G. 

[32] A.M. Ghez, M. Morris, E.E. Becklin, A. Tanner, and T. Kremenek, Nature, 407, 6802 (2000) 349-351. 

[33] M. Morris, Private communication. 

[34] A.M. Ghez et al., Astron.Nachr. 324 SI (2003) 527G. 

[35] P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schelke and E.A. Baltz, JCAP 0407 (2004) 008. 
[36] G. Gelmini, P. Gondolo A. Soldatenko C . Yaguna,|hep-ph/0605016 (2006) 
[37] G. Gelmini, P. Gondolo hep-ph/0602230 (2006) 

[38] L. Bergstrom, T. Bringmann, M. Eriksson, M. Gustafsson, Phys.Rev.Lett. 95, (2005) 241301. 

[39] S. Profumo, Phys.Rev.D 72, (2005) 103521. 

[40] J. F. Navarro, C. S. Frenk, S. D. M. White, ApJ 490 (1997) 493. 

[41] J. Silk, A. Stebbins, Astrophys.J. 411, (1993) 439. 

[42] E.W. Kolb, I. I. Tkachev, Phys.Rev.D 50, (1994) 769. 

[43] L. Bergstrom, J. Edsjo, P. Gondolo, and P. Ullio, Phys.Rev.D 59 (1999) 043506. 

[44] R. Aloisio, P. Blasi, A.V. Olinto, Astrophys.J. 601, (2004) 47. 

[45] A.V. Gurevich, K.P Zybin, Sov.Phys.Usp. 165, (1995) 723. 

[46] G. Bertone, D. Merritt, Mod.Phys.Lett.A 20, (2005) 1021. 

[47] M. Reid, ARA&A, 31, 345R (1993). 

[48] R. Genzel et al., ApJ 594, 812 (2003). 

[49] F. Aharonian et al., Nature 439 (2006) 695-698. 

[50] V.V. Vassiliev for the VERITAS Collaboration, Proc. 28th ICRC, Tskuba, Japan (2003). 
[51] L. Rolland and J. Hinton for the HESS Collaboration, Proc. 29th ICRC, Pune, India (2005). 



